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Quantum noise correlations have been employed in several areas in physics including condensed 
matter, quantum optics and ultracold atom to reveal non-classical states of the systems. So far, such 
analysis mostly focused on systems in equilibrium. In this paper, we show that quantum noise is also 
a useful tool to characterize and study the non-equilibrium dynamics of one dimensional system. 
We consider the Ramsey sequence of one dimensional, two-component bosons, and obtain simple, 
analytical expressions of time evolutions of the full distribution functions for this strongly-correlated, 
many-body system. The analysis can also be directly applied to the evolution of interference patterns 
between two one dimensional quasi-condensates created from a single condensate through splitting. 
Using the tools developed in this paper, we demonstrate that one dimensional dynamics in these 
systems exhibits the phenomenon known as " prethermalization" , where the observables of non- 
equilibrium, long-time transient states become indistinguishable from those of thermal equilibrium 
states. 

PACS: 67.85.-d 67.85. De 67.25.du 



I. INTRODUCTION 

Probabilistic character of Schrodinger wavefunctions 
manifests itself most directly in quantum noise. In many- 
body systems, shot-to-shot variations of experimental 
observables contain rich information about underlying 
quantum states. Measurements of quantum noise played 
crucial role in establishing nonclassical states of pho- 
tons in quantum optics pQ, demonstrating quantum cor- 
relations and entanglement in electron interferometers [2], 
and verifying fractional charge of quasi-particles in quan- 
tum Hall systems [3H5]. In atomic physics so far, noise 
experiments focused on systems in equilibrium. Re- 
cent work includes analysis of counting statistics in atom 
lasers[B], establishing Hanbury-Brown-Twiss effect for 
both bosons and fermionsjT], analysis of quantum states 
in optical lattices [8TU2j. observation of momentum corre- 
lations in Fermi gases with pairing [13] and investigation 
of thermal and quantum fluctuations in one and two di- 
mensional condensates [14II2T] . 

In this paper we demonstrate that analysis of quan- 
tum noise should also be a powerful tool for analyzing 
non-equilibrium dynamics of strongly correlated systems. 
Here we study the two equivalent dynamical phenomena; 
one given by the interaction induced decoherence dynam- 
ics in Ramsey type interferometer sequences for two com- 
ponent Bose mixtures in one dimension [55], and another 
given by the evolution of interference patterns of two one 
dimensional condensate created through the splitting of 
a single condensate [531 HI]- We obtain a complete time 
evolution of the full distribution function of the ampli- 
tude of Ramsey fringes or interference patterns. In the 
case of Ramsey fringes, the average amplitude of Ram- 
sey fringes measures only the average value of the trans- 



verse spin component. On the other hand, full distribu- 
tion functions are determined by higher order correlation 
functions of the spins. Hence full distribution functions 
contain considerably more information about the time 
evolution of the svstem [24H27] and provide a powerful 
probe for the nature of the quantum dynamics under 
study. In particular, we use the simple expressions of 
full distribution functions to demonstrate the phenomena 
of "prethermalization" in these one dimensional systems, 
where observables in non-equilibrium long-time transient 
states become indistinguishable from those in thermal 
equilibrium states. 

One dimensional systems with continuous symmetries, 
including superfluids and magnetic systems, have a spe- 
cial place in the family of strongly correlated systems. 
Quantum and thermal fluctuations are so extreme that 
long range order is not possible in equilibrium. Such 
systems can not be analyzed using standard mean-field 
approaches, yet they can be studied through the applica- 
tion of methods specific to one dimension such as exact 
Bethe ansatz solutions [28 39], effective description us- 
ing Tomonaga-Luttinger and sine-Gordon models [3DH1S], 
and numerical analysis using density-matrix renormal- 
ization group(DMRG) and matrix product state (MPS) 
methods [46]. Such systems are often considered as gen- 
eral paradigms for understanding strongly correlated sys- 
tems. One dimensional systems also give rich examples of 
integrable systems, where due to the existence of infinite 
number of conserved quantities, equilibration does not 
take place |47tf4"9"] . Hence the problem we consider in this 
paper is important for understanding fundamental issues 
such as the quantum dynamics of strongly correlated sys- 
tems and equilibration/non-equilibration of many-body 
systems, as well as for possible applications of spinor 



condensates in spectroscopy, interferometry, and quan- 
tum information processing [501452] . 



Our work is motivated by recent experiments of 
Widera et al.^?Z\ who used two hyperfine states of 87 Rb 
atoms confined in 2D arrays of one dimensional tubes 
to perform Ramsey type interferometer sequences. They 
observed rapid decoherence of Ramsey fringes and the 
near absence of spin echo. Their results could not be 
explained within the single mode approximation which 
assumes a macroscopic Bose condensation into a single 
orbital state, but could be understood in terms of the 
multi-mode Tomonaga-Luttinger type model. Yet the 
enhanced decoherence rate and suppression of spin echo 
do not provide unambiguous evidences for what the ori- 
gin of decoherence is. In this paper, we suggest that the 
crucial evidence of the multi-mode dynamics as a source 
of decoherence should come from the time evolution of 
the full distribution functions of the Ramsey fringe am- 
plitude. Such distribution functions should be accessible 
in experiments on Atom Chips 14, 53 56 from the anal- 
ysis of shot-to-shot fluctuations. 

This paper is organized as follows. In Section [TTJ 
we describe two physically distinct, but yet mathemat- 
ically equivalent, dynamics in one dimensional systems, 
namely, the dynamics of spins in a Ramsey sequence and 
the dynamics of phase and contrast in interference pat- 
terns between two split condensates created from a single 
condensate. We start with illustrating the basic physics 
governing the dynamics studied in this paper, and give a 
summary of the central results including the prediction 
of prethermalization phenomena. The formal descrip- 
tions of the details of the theory for spin dynamics in 
Ramsey sequence start from Section III where we give 



the Hamiltonian of the one dimensional system based on 
Tomonaga-Luttinger approach. In Section [iVj we derive 
the analytical expression for the time evolution of the full 
distribution function for a simple case in which charge 
and spin degrees of freedom decouple. This decoupling 
limit gives a good approximation to the experimental sit- 
uation of Widera et al\22\ . A short summary of the result 
in this decoupling limit has been already reported in Ref 
[52j . More general case in which spin and charge degrees 
of freedom mix is studied in Section [V] Such mixing in- 
troduces the dependence of spin distribution functions on 
the initial temperature of the system. All the results ob- 
tained in Sections [IV| and |V| can be extended to the study 
of the dynamics of interference patterns, using the map- 
ping described in Section [TTJ The details of the dynamics 
of phases and contrasts in interference patterns between 
split condensates is studied in Section |VI| We demon- 
strate the prethermalization phenomena and show that 
the interference contrasts of split condensates in a steady 
state have indistinguishable distributions from those of 
thermal condensates at some effective temperature T e jf. 
We conclude in Section [Vll| with a discussion of possible 
extensions of this work. 
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FIG. 1: Ramsey sequence for one dimensional system with 
two component bosons considered in this paper. (1) All atoms 
are prepared in spin up state; (2) 7r/2 pulse is applied to rotate 
each atom into the x direction; (3) spins freely evolve for time 
t. In actual experiments, final 7r/2 pulse is applied to measure 
the x component of spin. The imaging step (4) is omitted in 
the illustration. In this paper, spin operators refer to the ones 
before the final n/2 pulse. 



II. DESCRIPTION OF ONE DIMENSIONAL 
DYNAMICS AND THE SUMMARY OF RESULTS 

A. Ramsey Dynamics 

In this paper, we study the dynamics of one dimen- 
sional, interacting two-component Bose mixtures in the 
Ramsey-type sequence. In analogy with spin-1/2 parti- 
cles, we refer to one component to be spin-up and the 
other component to be spin-down. In the experiment 
of cold atoms in Ref [22] , two hyperfine states are used 
for these two components. In the following, we consider a 
generic situation where there is no symmetry that relates 
spin-up and spin-down. In particular, unlike fermions 
with spin-1/2, there is no SU(2) symmetry. In a typical 
experimental setup with cold atoms, there is a harmonic 
confinement potential along the longitudinal direction of 
condensates, but here we assume the absence of such a 
harmonic trap potential. Our consideration gives a good 
approximation for the central region of cold atom exper- 
iments in the presence of such potentials. 

Ramsey- type sequence is described as follows (Figure 
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1. All atoms are prepared in spin up state at low tem- 
perature 

2. 7r/2 pulse is applied to rotate the spin of each atom 
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into the x direction 

3. Spins evolve for time t 

4. Spins in the transverse direction (x — y plane) are 
measured 

In a typical experimental situations |22j. the last mea- 
surement step is done by applying a 7r/2 pulse to map 
the transverse spin component into z direction, which 
then can be measured. In the following discussions, we 
describe the dynamics in the rotating frame of Larmor 
frequency in which the chemical potentials of spin-up 
and spin-down are the same in the absence of interac- 
tions. In this frame, the evolution of spins in the third 
step is dictated by the diffusion dynamics coming from 
interactions. Unlike the conventional use of the Ramsey 
sequence in the context of precision measurements, here 
we employ the Ramsey sequence as a probe of correlation 
functions in one dimensional system. 

The description of the spin dynamics starts from the 
highly excited state prepared after the tt/2 pulse of step 2. 
The subsequent dynamics during step 3 crucially depends 
on the nature of excitations in the system. In particular, 
the dynamics of two-component Bose mixture in one di- 
mension is quite different from those in three dimension. 
In three dimensions, bosons form a Bose-Einstcin con- 
densate(BEC) at low temperature, and particles occupy 
a macroscopic number of k = mode. Then, the spin 
diffusion of three dimensional BEC is dominated by the 
spatially homogeneous dynamics coming from the single 
k = mode at sufficiently low temperatures. On the 
other hand, bosonic systems in one dimension do not 
have the macroscopic occupancy of the k = mode, and 
their physics is dominated by the strong fluctuations, to 
the extent that the system cannot retain the long range 
phase coherence even at zero temperature [H]. Thus, the 
spin dynamics of one dimensional bosonic system nec- 
essarily involves a large number of modes with different 
momenta and the spin becomes spatially inhomogeneous 
during the step 3 above. 

Such dynamics unique to one dimension can be probed 
through the observation of transverse spin components in 
the fourth step. Since we aim to capture the multi-mode 
nature of the dynamics in one dimension, we consider the 
observation of spins at length scale I, given by 
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drS a (r,t) 
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where S°(r, t) with a = x,y are the transverse compo- 
nents of spin operators after time evolution of step 3 of 
duration t. We assume that I is much larger than the spin 
healing length £ s , and much smaller than the system size 
L to avoid finite size effects. Furthermore, we assume 
that the number of particles within the length I, Ni, is 
large, so that the simultaneous measurements of Sf(i) 
and Sf(t) are in principle possible. For large Ni, the 
non-commutativity of Sf and Sf gives corrections of the 



order of 1/y/Ni compared to the average values. In this 
situation, it is also possible to measure the magnitude 

of transverse spin components, S^ = <J {Sf} + (jSf 

which we will extensively study in the later sections. 

Due to quantum and thermal fluctuations, the mea- 
surements of Sf (i) give different values from shot-to-shot. 
After the 7r/2 pulse of step 2, the spins are prepared in x 

direction, so the average value yields (Sf(t = Q)\ ~ Ni/2 

and (sf(t = 0)^> ~ 0. In the rotating frame of Larmor 
frequency, the subsequent evolution does not change the 
expectation value of the y component so that ( Sf (t) 



throughout. The decay of the average ^Sf(t)J during 
the evolution in step 3 tells us the strength of spin dif- 
fusion in the system. The behaviors of /Sf due to 

spin diffusion are similar for one and three dimensions, 
and the difference is quantitative, rather than qualita- 
tive. On the other hand, a richer information about the 
dynamics of one dimensional system is contained in the 
noise of Sf (t) . Such noise inherent to quantum systems is 
captured by higher moments ( ( Sf(t) )"). In this paper, 
we obtain the expression for the full distribution func- 
tion P ; a (a,i) which can produce any moments of Sf(t) 
through the relation 



<(SfW 



= J dQP,"(Q,t)a n , 



(2) 



where P"(a, t)da represents the probability that the 
measurement of Sf(i) gives the value between a and 
a + da. We will see in Section [TVT t hat it is also possible 
to obtain the joint distributions P l ' y (a, (3, t) of Sf (t) and 
Sf(t) as well as the distribution P [ ± (a,t) of the squared 

transverse magnitude ■ 

Now we summarize the main results of this paper, and 
give a qualitative description of spin dynamics in the 
Ramsey sequence. Elementary excitations of spin modes 
in the system are described in terms of linearly dispers- 
ing spin waves with momenta k and excitation energies 
c s \k\, where c s is the spin wave velocity. When certain 
symmetry conditions are satisfied(see discussion in Sec- 
spin and charge degrees of freedom decouple, 



IV) 



tion 

and these spin waves are free and they do not interact 
among themselves in the low energy descriptions within 
so-called Tomonaga-Luttinger theory [41] 02]. Here, we 
describe the result in this decoupling limit, but the qual- 
itative picture does not change even after the coupling 
between spin and charge is introduced, as we will see in 
Sec. El 

The initial state prepared after ir/2 pulse in step 2 in 
which all spins point in the x direction is far from the 
equilibrium state of the system because interactions of 
spins are not symmetric in terms of spin rotations. Thus, 
the initial state contains many excitations and the subse- 
quent dynamics of spins is determined by time evolution 
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FIG. 2: Illustration of the dynamics of spins in the presence 
of spin wave excitations. At short times(top), high momenta 
excitations contribute to fluctuations of the spins, but their 
effect is weak. At long times(bottom), low momenta excita- 
tions lead to the strong fluctuations of the spins. Such fluc- 
tuations with wavelengths larger than I rotate the regions of 
length I as a whole so that they do not lead to the decay of 
the magnitude of spin §i , but result in diffusion of Sf . 



of the spin waves. A spin wave excitation with momen- 
tum k rotates spins with length scale ~ 2ir/k and time 
scale ~ l/(c s |fc|). The amplitude of fluctuations coming 
from the spin wave with momentum k is determined by 
the initial state as well as the nature of spin wave ex- 
citations. We find that the energy stored in each mode 
is approximately the same(see discussion in Sec. IV B 5 1, 



thus the amplitude of fluctuations for wave vector k scales 
as 1/k 2 . Therefore, the fluctuation of spins is weak at 
short wave lengths and short times, and strong at long 
wave lengths and long times. In Fig(2] we illustrate such 
dynamics of spins due to fluctuations of spin wave exci- 
tations. It leads to the distributions presented in Figj3] 
and Figj4j Here, we have plotted distribution function of 

the squared transverse magnitude of spins {sf^j (Fig 3) 

and the joint distribution function(Figj4]) with L = 200, 
K s = 20 and various integration length l/£ a = 20, 30, 40. 
K s is the spin Luttinger parameter, which measures the 
strength of correlations in ID system (see Eq. Q below) , 
and £ s is a spin healing length which gives a characteristic 
length scale in the low energy theory of spin physics. 

The multi-mode nature of one dimensional system, in 
which spin correlations at different length scales are de- 
stroyed in qualitatively different fashion during the dy- 
namics, can be revealed most clearly in the squared trans- 
verse magnitude of spins (jS^ , plotted in Fig 
initial state, all the spins are aligned in the x direction, 
so the distribution of (•S/ 1 ) i s a delta function peak 

at its maximum value, ~ (pi) 2 , where p is the average 
density of spin-up or spin-down. The evolution of spin 
waves lead to the fluctuations of spins and thus to the 
decay of the integrated magnitude of the transverse spin. 
How the spin waves affect the integrated magnitude of 
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FIG. 3: Time evolution of the distribution P l L (a) of the 
squared transverse magnitude of spins, ^S^) , f° r the system 
size L/^ a = 200, the spin Luttinger parameter K s = 20 and 
various integration length Z/£ s = 20,30,40. Here £ s is the 
spin healing length, and the x axis is scaled such that the 
maximum value of a is 1. Time is measured in units of £ s /c s 
where c s is the spin sound wave velocity. The evolution of 
the distribution crucially depends on the integration length. 
The steady state of the distribution of the squared transverse 
magnitude has a peak at a finite value for short integration 
length l/^s = 20, whereas the peak is at for long integration 
length I /l s = 40. 



spins strongly depends on the wavelength of the excita- 
tions. Spin excitations with momenta much smaller than 

~ 2tt/Z do not affect smce these spin waves ro- 

tate the spins within I as a whole, while spin excitations 
with higher momenta lead to the decay of the magnitude. 
This is in stark contrast with the x component of the spin 
Sf, which receives contributions from spin waves of all 
wavelengths. 

As a result of different contributions of spin wave exci- 
tations with different wavelengths to the integrated spin 
magnitude, there are two distinct behaviors of the dis- 
tributions of > onc f° r short integration length I, 
which we call " spin diffusion regime" and another for long 
integration length I, which we call "spin decay regime." 

For short integration length I, the distribution func- 
tion of i s always peaked near its maximum value 

(pi) 2 during the dynamics because the strengths of fluc- 
tuations coming from spin waves with high momenta are 
suppressed by l/fc 2 (Fig[|} l/£ a = 20). While the magni- 
tude of spins does not decay in this regime, fluctuations 
still lead to a diffusion of Sf, thus, we call this regime 
the "spin diffusion regime." 

On the other hand, for long integration length, spin 
waves lead to fluctuations of the spins within the inte- 
gration region, and the spins are randomized after a long 
time. This randomization of spins leads to the devel- 
opment of a Gaussian-like peak near (S^) 2 — 0(Figji 
l/£ s = 40). During the intermediate time, both peaks a ; 
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FIG. 4: Time evolution of the joint distribution function P x ' y (a,/3) for the system size L/£ s = 200, the spin Luttinger 
parameter K s = 20 and various integration lengths l/£ s — 20(left),30(middle), and 40(right). Time is measured in units of 
£ s /c s where c s is the spin sound wave velocity. Here axes are scaled such that the maximum value of a and /3 are 1. For for short 
integration length l/£ 3 — 20, the dynamics leads to the distribution with the "ring"-like structure, showing that the magnitude 
of spins does not decay much (spin diffusion regime). On the other hand, for longer integration lengths, the magnitude of spins 
decays quickly and the distribution forms a "disk"-like structure(spin decay regime). 



and the maximum value (pi) 2 are present, and one can 
observe the double peak structure. Because of the strong 
decaying behavior of the magnitude of spin, we call this 
regime the "spin decay regime." 

More complete behaviors of distribution functions can 
be captured by looking at the joint distribution func- 
tions from ^which we can read off the distributions of 

both (§fj and Sf , see Fig 4J57]. In the "spin diffusion 

regime" with short I, the joint distributions form a "ring" 
during the time evolution, whereas in the "spin decay 
regime" with long I, they form a "disk" -like structure in 
the long time limit. As we will see later, a dimensionless 

2 i 

parameter given by ~ 4 ^- 1 determines whether the 
dynamics belongs to the "spin diffusion regime" (lo < 1) 
or the "spin decay regime" (lo ^ 1). 

We emphasize that in three dimensions, spin waves 
are dominated by k = mode and therefore, there is al- 
most no decay in the magnitude of spins throughout the 
dynamics. Therefore, the existence of two qualitatively 

different behaviors of distribution functions of f 1 ^/ - ) un ~ 
ambiguously distinguishes the dynamics in one and three 
dimensions. 
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FIG. 5: The interference of two quasi-condensates that are 
created by splitting one single quasi-condensate. After the 
splitting, the quasi-condensates are held for time t and then 
the transverse confinement is released. The two condensates 
interfere with each other after the release and the position of 
constructive interference is denoted by solid line in the figure. 
This interference pattern contains the information about the 
local phase difference between the two quasi-condensate at 
the time of release. 



B. Dynamics of interference between split 
condensates 

The dynamics of Ramsey sequence considered above 
can be directly mapped to the dynamics of interference 
pattern of split one dimensional quasi-condensate |58j . 
More specifically we consider the following sequence of 
operations(see Figj5]). First, we prepare one-component 
ID quasi-condensate in equilibrium. At time t = 0, the 
quasi-condensate is quickly split along the axial direction, 



and the resulting two quasi-condensates are completely 
separated. The two quasi-condensates freely evolve for 
a hold time of t, and they are released from transverse 
traps to observe the interference pattern between these 
two quasi-condensates. Such dynamics of the interference 
patterns as a function of hold time t has been observed in 
the experiment by Hofferberth et. al. |23j . where the aver- 
age of interference patterns is analyzed in details. Here 
we study the unique one dimensional dynamics by look- 
ing at the full distributions of interference patterns (for 
an experimental study see: Gring et al|24). ). 
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The dynamics of split condensates can be mapped to 
the dynamics of the Ramsey interferometer studied in 
this paper. The splitting of a quasi-condensate corre- 
sponds to the initial 7r/2 pulse in the Ramsey sequence. 
If we call one of the quasi-condensates to be L for left 
and another to be R for right, L(i?)-condensate corre- 
sponds to spin-up (spin-down) component. Thus, the 
density difference between the two condensates corre- 
sponds to z component of the spin, namely S z (r,t) = 
ip^ L {r, t)ij)i,{r, t) — ily R {r, t)tpn{r, t) where ipl(r, t) is the cre- 
ation operator of particles in i = L, R(R) condensate. 
Moreover, the local phase difference between the two con- 
densates corresponds to the local spin direction in x — y 
plane. To see this, we first note that ?/>|(r, t)ipn(r, t) cor- 
responds to the spin raising operator S + (r, t). This oper- 
ator is expressed in terms of the phase difference between 
the condensates 4> s (r,t) as S + (r, t) = ^ L {r,t)Tj)R(r,t) ~ 
p e i<f> s (r,t) w ] iere p [ s the average density of each conden- 
sate. Thus, for example, x and y spin operators are given 
by S x (r,t) = pcos(0 s (r,i)) and S»(r,t) = psm($ s (r,t)). 
Immediately after the splitting, the phases of the two 
quasi-condensates at the same coordinate along the axial 
direction are the same. Therefore, the splitting prepares 
spins in x direction in the language of the Ramsey se- 
quence, and thus the splitting effectively amounts to the 
7r/2 pulse. 



The interference of two quasi-condensates measures the 
local phase difference at time t. If the phases of L and R 
condensates are the same, the interference pattern has 
a constructive peak at the center of two condensates, 
which we call x — 0. Thus, a shift in the interfer- 
ence pattern measures the local phase difference between 
the two condensates, which yields the information about 
S x (r) = pcos(cj) s (r)) as well as S v (r) = psm(cj) s (r)). We 
note that the integrated interference contrast that can 
be extracted from experiments is given by the expres- 

2 



sion C 2 = 



Pfl 



i(p s (r) 



and related to the transverse 



magnitude of spins as (jfij 



Sf + iSf 



= C 2 



All the results obtained for Ramsey dynamics are di- 
rectly applicable to the dynamics of interference patterns 
between split condensates. When the splitting process 
prepares two condensates with equal average number of 
particles, "spin" and "charge" degrees of freedom de- 
couple, see Sec |VI| for details. In particular, depend- 
ing on the integration length of the interference patterns 
along the axial direction, there exists two regimes corre- 
sponding to "phase diffusion regime" and "contrast decay 
regime," analogous to "spin diffusion regime" and "spin 
decay regime" described in the previous section, respec- 
tively. 



C. Prethermalization of one-dimensional 
condensates 

The equilibration and relaxation dynamics of generic 
many-body systems are fundamental open problems. 
Among possible processes, it has been suggested that 
the time evolution of some systems prepared in non- 
equilibrium states results in the reaching of quasi-steady 
states within much shorter time than equilibration time. 
This quasi-steady state is often not a true equilibrium 
state, but rather it is a dephased state, and true equi- 
libration takes place at much longer time scale. Yet in 
some cases, the physical observables in the quasi-steady 
states take the value corresponding to the one in thermal 
equilibrium at some effective temperature T e g, display- 
ing disguised "thermalized" states. Such surprising non- 
equilibrium phenomena, called prethermalization, have 
been predicted to occur in quantum as well as classi- 
cal many-body systems [59- BT] and observed in Gring et 

aim- 

In particular, integrable one dimensional systems are 
known not to thermalize and indeed, experiments in 
Ref. [UJ have observed an exceedingly long equilibra- 
tion time. Yet even in this extreme case, we suggest in 
this section that many-body one dimensional systems can 
reach disguised "thermalized" states through prethermal- 
ization phenomena within a short time. 

Ramsey dynamics and dynamics of interference pat- 
terns between split condensates described in the previous 
sections are particular examples of dynamics in which 
slow equilibration is expected because the system essen- 
tially consists of uncoupled harmonic oscillators in the 
low energy description (see Sec IV I . In the following, we 
give a heuristic argument that the distribution of the in- 
terference contrast amplitudes of the two non- equilibrium 
quasi-condensates are given by that of two equilibrium 
quasi-condensates at some effective temperature T e g. We 
give more details in Sec. |VIC| 

Long time after the splitting, the position of the inter- 
ference peak becomes completely random, and yields no 
information. Therefore, we focus on the squared trans- 
verse magnitude of the spin ("Sf^ > or equivalently, the 

interference contrasts C 2 . In the following, we describe 
the physics for the dynamics of split condensates, but 
the same argument can be applied to Ramsey dynamics. 
The interference contrast of the split condensates after a 
long time is determined by the "average" phase fluctua- 
tions present in the system. In the dephased limit, such 
fluctuations are determined by the total energy present 
in each mode labeled by momenta k. Now for sufficiently 
fast splitting, the energy -E S piit contained in each mode is 
independent of momenta because the density difference 
of quasi-condensates along the axial direction is uncorre- 
cted in the initial state beyond the spin healing length 
£ s (see discussion below Eq. (10)). On the other hand, 
the interference contrast of thermal condensates at tem- 
perature T is determined by the thermal phase fluctua- 
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tions caused by excitations whose energy is distributed 
according to equipartition theorem; each mode in the 
thermal condensates contains the equal energy of kgT. 
Thus from this argument, we find that the interference 
contrast of split condensates after a long time becomes in- 
distinguishable from the one resulting from thermal con- 
densates at temperature ^T e g = E sp n t . We will show 
in Sec | VI C| that in fact, the full distribution function of 
interference contrast becomes indistinguishable for these 
two states. In the case of splitting two dimensional con- 
densates, equipartition of energy and existence of "non- 
equilibrium temperature" was pointed out in Ref. |61|. 

Here we propose the occurrence of such prethermaliza- 
tion within Tomonaga-Luttinger theory. We emphasize 
that within Tomonaga-Luttinger theory of low energy 
excitations, different modes are decoupled and there- 
fore no true thermalization can take place. In realistic 
experimental situations, such integrability can be bro- 
ken and relaxation and thermalization process are ex- 
pected to occur after a long time dynamics. The re- 
quirement to observe the prethermalization phenomena 
predicted in this theory in experiments depends on the 
time-scale of other possible thermalization processes we 
did not consider in our model such as the effective three- 
body collisions [BH [S3], relaxation of high energy quasi- 
particles[64], or interactions among the collective modes 
through anharmonic terms we neglected in Tomonaga- 
Luttinger theory [55], When all these processes occur at 
much slower time scale than the prethermalization time- 
scale, which is roughly given by the integration length 
divided by the spin sound velocity ~ l/c s for decoupling 
case, the observation of prethermalization should be pos- 
sible. In one dimension, the dynamics is strongly con- 
strained due to the conservation of energy and momen- 
tum, and therefore it is likely that the dynamics is dom- 
inated by the modes described by Tomonaga-Luttinger 
theory for long time for quasi-condensates with low ini- 
tial temperatures. 

Such prcthcrmalizations are expected to occur even in 
higher dimensional systems [TO1 [25, 61 ( 66 . We note that 
the conditions for the experimental observations of the 
phenomena might be more stringent because true ther- 
malization processes are expected to take place much 
more quickly in two and three dimensional systems. 



III. TWO COMPONENT BOSE MIXTURES IN 
ONE DIMENSION: HAMILTONIAN 



In this paper, we study the dynamics of two-component 
Bose mixtures in one dimension through Tomonaga- 
Luttinger formalism [41] [42]. As we have stated before, 
we assume the rotating frame, in which spin-up and spin- 
down particles have the same chemical potential in the 
absence of interactions. The Hamiltonian of two compo- 



nent Bose mixtures in one dimension is given by 

i-L/2 



dr V]— V$(r)Wi(r) 
-l/2 * 2m l 



(3) 



Here ipi with i =t,| describe two atomic species with 
masses m, and gij are the interaction strengths given by 
9ij = v ' a *j [57] where v± is the frequency of transverse 
confinement potential and are the scattering lengths 
between spin i and j. System size is taken to be L, and 
we take the periodic boundary condition throughout the 
paper. In addition, we use the units in which h = 1. 

In the low energy description, the Hamiltonians for 
weakly interacting bosons after the initial tt/2 rotation 
can be written in quadratic form, and given by 



H 


— Hf + H_i 




,L/2 




= / dr 


J -L/2 




l-L/2 


Hi 


= dr 


J -L/2 




fL/2 


Hint 


= 21 a 




J -L/2 



(4) 



2m^ 



2m; 



(V</» T (r)) 2 +.g n (n t (r)) 2 



where p is the average density of each species and n a 
are variables representing the phase and density fluc- 
tuation for the particle with spin a. These variables 
obey a canonical commutation relation [n a (r), <p a (r')] — 
—i5(r — r'). In the Hamiltonian above, we included the 
kinetic interaction term g^, which is zero for weakly in- 
teracting bosons, but allowed by inversion symmetry and 
non-zero for generic Tomonaga-Luttinger Hamiltonians. 

We note that in the weakly interacting case, one can 
obtain the parameters of the Hamiltonian in Eq. Q 
such as and nii through hydrodynamic linearization 
of the microscopic Hamiltonian. In this case, we assume 
the small fluctuations of the densities h a and phases 
V0 a and expand the Hamiltonian in Eq. {[3J) to the 
second order in these variables through the expression 
ijp a ~ y/p + h a e % ^" , resulting in the form of the Hamilto- 
nian in Eq. Q . Due to this assumptions of small spatial 
variations of the phase, h a and <f> a represent the "coarse- 
grained" variables where collective modes have linear 
dispersions. The Hamiltonian of Tomonaga-Luttinger 
theory in Eq. Q can also describe effective low-energy 
physics of strongly interacting systems, but in this case 
there is no simple relation between microscopic parame- 
ters and the parameters of the Hamiltonian in Eq. Q. 

In order to describe the spin dynamics, we define spin 
and charge operators as the difference and the sum of spin 
up and down operators, i.e. 4> s — 4>^ — 4>±, <f> c = <f>^ + <j>±, 
n s = \{n\— n c = ^(h^ + fii). In this representation, 
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the Hamiltonian in Eq. Q becomes 



H — H s + H c + H m i X 

rL/2 

H s = / dr 

J-L/2 
r L/2 

If,. I dr 

-L/2 
rL/2 



P 

2m, 

P 

2m, 



(V0 s (r)) 2 +.g s (n s (r)) 2 
(WMr)) 2 +g c (h c (r)) 2 



(5) 
(6) 



dr 



L/2 



g m ixn s (r)n c (r) + .g^WsMWcM 

(7) 



where interaction strengths are given by g c = g^ + g^ + 

2 514' 9s = 9ff + 9 '44 - 2,g^, g mix = g^ - g±±, g mix = 
pj (8to^) — p/(8m±). The masses are given by the relations 
p/(2m c ) = p/(8m t ) + p/(8m ; ) + gfj2 and p/(2m s ) = 
p/(8m T )+p/(8m ; )-^/2. 

The spin variables <p s and n s are "coarse-grained" in 
the sense that they represent the operators in the long 
wavelength beyond the spin healing length £ s . £ s is deter- 
mined from microscopic physics and gives the length be- 
low which the kinetic energy of spins wins over the inter- 
action energy, see Eq. ([5| above. For weakly interacting 
bosons with m-j- = m\, = m, it is given by £ s = it / ^Jmpg s . 
In the following, we assume that the number of particles 
within the spin healing length is large, i.e. £ s p ^> 1. 
This condition is always satisfied for weakly interacting 
bosons. 

In the next section, we consider the case g m ix = and 
9 mix = 0, in which spin and charge degree of freedom 
decouple. Then the dynamics of spins is completely de- 
scribed by the spin Hamiltonian in Eq. The general 
case in which g m i X ^ and g mix ^ will be treated in 
SeclVl 



IV. DYNAMICS OF FULL DISTRIBUTION 
FUNCTION FOR DECOUPED SPIN AND 
CHARGE DEGREES OF FREEDOM 

A. Hamiltonian and initial state 

The experiment of Widera et al. [25] used F — 1, mp = 
+1 and F = 2,mp = —1 states of 87 Rb for spin-up and 
spin-down particles, respectively. These hyperfine states 
have the scattering lengths a aa * such that a++ ~ a^. 
Consequently, the mixing Hamiltonian in Eq. ufh approx- 
imately vanishes for weak interactions. Motivated by this 
experiment, here we consider the decoupling of spin and 
charge degrees of freedom [52]. Spin dynamics in this case 
is completely determined by the spin Hamiltonian 



H s = - 
2 



A", 



-(V<Mr)) 2 



n 2 s (r) 



K, 



dr (8) 



where K s is the spin Luttinger parameter representing 
the strength of interactions, and c s is spin sound velocity. 



K s and c s are directly related to the spin healing length 
£ s in the weak interaction limit, given by 2K S = p£ s 
and c s = 2 m £ ■ ™ s ( r > ^) * s ^ ne l° ca l spin imbalance 
h s = ^Ui^afl)^ an< ^ 4>s{ r it) is related to the direc- 
tion of the transverse spin component pe l ^ s = 4''L 'ap' t l ; P- 
Here, -01 is the creation operator of spin a =f , 4. These 
variables n s and 4> s obey a canonical commutation rela- 
tion [h s (r), 4> s (r')] = —iS(r -r'). 

Other spin variables can be similarly defined in terms 
of coarse grained spin variables h s and 4> s . In the follow- 
ing, we consider the general transverse spins pointing in 
the direction (x,y,z) — (cos 9, sin 9,0) integrated over I 
given by 



St = 



1/2 



dri/jl(r,t) cos. 

-1/2 \ 2 

1/2 



dr- 



1/2 



-i(4> s (r)-6) 



(9) 



where <r a with a = x,y are Pauli matrices. Here Sf 
with 9 = corresponds to spin x operator and 9 = n/2 
corresponds to spin y operator. In order to explore the 
one dimensional dynamics resulting from Hamiltonian in 
Eq.([8]), we analytically compute the mth moment of the 

spin operator Sf, ((§1^ ^; after time t of the ir/2 pulse 
of the Ramsey sequence. Then, the full distribution func- 
tions of Sf and Sf, as well as the joint distribution of 

these will be obtained from ( ( S, 



In order to study the dynamics of Ramsey interferom- 
eter in terms of low energy variables n s and <f> 3 , we need 
to write down an appropriate state after the it/2 pulse 
in terms of h s and <j> s . If pulse is sufficiently strong, 
each spin is independently rotated into x direction af- 
ter the 7r/2 pulse. Naively, this prepares the initial state 
in the eigenstate of S x (r) = p cos <j> s {r) with eigenvalue 
<j) s {r) = 0. However, due to the commutation relation 
between n s and <j> s , such an initial state has an infinite 
fluctuation in h s and therefore, the state has an infi- 
nite energy according to Eq.Q. This unphysical con- 
sequence comes about because the low energy theory in 
Eq. ^ should not be applied to the physics of short time 
scale given by 1 / E c where E c is the high energy cutoff of 
Tomonaga-Luttinger theory. During this short time dy- 
namics, the initial state establishes the correlation at the 
length scale of spin healing length £ s . The state after this 
short time dynamics can now be described in terms of the 
coarse-grained variables n s (r) and 4> s (r). The variables 
n s (r) and 4> s (r) are defined on the length scale larger 
than the spin healing length £ s . Since the z component 
of spins are still uncorrelated beyond £ s after the initial 
short time dynamics, the appropriate initial condition of 
the state is written as 



(S z (r)S z (r')) - (n s (r)n s (r')) 



PV 



S(r 



(10) 
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where the delta function 6(r — r') should be understood 
as a smeared delta function over the scale of £ s . Because 
the state after the short time dynamics is still close to the 
eigenstate of the S x (r) operator, spins are equal superpo- 
sitions of spin-up and spin-down. Then the distribution 
of Sf — / S z {r)dr is determined through random pick- 
ing of the values ±1/2 for 2pl particles. Due to the cen- 
tral limit theorem, the distribution of Sf = S z (r)dr is 

Gaussian, i.e. ( ( Sf \ ) — (plrj / 2) n . In particu- 



lar, / (s?J ) = which determines the magnitude 



of the fluctuation for S z (r) in Eq. ( 10 ) . In Eq. ( 10 ) , we also 



introduced the phenomenological parameter i] which ac- 
counts for the decrease and increase of fluctuations com- 
ing from, for example, imperfections of ir/2 pulse. The 
ideal, fast application of tt/2 pulse corresponds to r/ = 1. 
In the experimental realization of Ref. |22| , r\ was deter- 
mined to be between 0.8 and 1.3 through the fitting of 
the experiment with Tomonaga-Luttinger theory for the 
time evolution of the average x component of the spin, 
(j3f\. Through engineering of the initial state such as 

the application of a weak tt/2 pulse, rj can also be made 
intentionally smaller than 1. 

A convenient basis to describe the initial state of the 
dynamics above is the basis that diagonalizes the spin 
Hamiltonian of Eq.dsl). The phase and density of the 
spins 4> a (r) and n s (r) can be written in terms of the cre- 
ation b k and annihilation b Si k operators of elementary 
excitations for the spin Hamiltonian in Eq.d8|) as 



by 



\fc#0 / 



(n a ,o|V'«,fc=o) = exp 



2pri 



(12) 



where 2Wk 



1 — Qfc 



\k\K a 



Here the state \n s 



l+afc ' *■* -Kp-q 

is the normalized eigenstate of the operator n Sj o with 
eigenvalue n Sj o- The summation of k in the exponent 
has a ultraviolet cutoff around k c — 2tt/^ s . AT is the 
overall normalization of the state. It is easy to check 
that (0o | n Sl fen Sl fe' \ipo) = %>$k,-k', which corresponds to 



Eq. (10). 



B. Moments and full distribution functions of spins 

After free evolution of the initial state \ipo) for time 
t, the state becomes \tp(t)} = e~ %Hst \tpo)- We charac- 
terize the state at time t by the mth moments of spin 

operators, ^(Sf ) m ^- As we will see below, the full distri- 
bution function can be constructed from the expression 



of as, 



?8\m 



181 



We consider the evaluation of moments ( (S"Y 



at 



time t, \ip(t)}. Each momentum k component of the ini- 
tial state \ipo) independently evolves in time. Since k = 
mode has a distinct evolution from other k ^ modes, 
we separately consider k — and k ^ modes. 



Vl 



//,-/• 



(b] k - b s ^ k y kr + S , O 



h s (r) 



It ^ 



Vl 



n s he 



ikr 



J2c s \k\bl k b S:k + 

fc#0 



-8 -2 



2K, 



l s,0 



(11) 



where we defined 



and n s f, to be the Fourier trans- 



form of operators 4> s (r) and h s (r) . 6^ . creates a col- 
lective mode with momentum k and follows a canonical 
commutation relation [b Sy k,bl k ] — 1. Note that k = 
mode has no kinetic energy, and it naturally has differ- 
ent evolution from k ^= modes. 



The Gaussian state determined by Eq.(10) takes the 



form of a squeezed state of operators b Si k, and it is given 



1. k — mode 



The Hamiltoni an o f k — mode is given by H. 



2K B "«,o 



in Eq. jTTj ). Therefore, in the basis of 



,fc=0 



k = part of the state \ip(t)} is given by 



{n s ,o\ e 



-iH, 



{(- 



k = 0t \4>k=0 

L • 7TC S t 



' 2K S 



) n 8,0}. 



(13) 



where A4=o is the normalization of the state. The initial 
Gaussian state of n S) o stays Gaussian at all times, and 
any analytic operator of (j3s,o and n Sr o can be exactly 
evaluated through Wick's theorem. For example, k = 
part contributes to the decay of the average of the x 



component of spin ( S 



fc=0 



sr 



— Ipe 
1 



lpRe^( 



c s irt 

~k7 



VP 
2 



(14) 



This diffusion of the spin from k = contribution is gen- 
erally present in any dimensional systems, and not par- 
ticular to one dimension. Physical origin of this diffusion 
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is the interaction dependent on the total spin , S% , The 
eigenstate of S x with eigenvalue pi is the superposition of 
different eigenstates of S z with eigenvalues m z , and they 
accumulate different phases e~ m * in time. This leads to 
the decay of (^S^j. In the thermodynamic limit L — > oo, 
the uncertainty of S z becomes diminishingly small, and 
therefore, the decay of coming from k = goes 

to zero. More interesting physics peculiar to one dimen- 
sional systems comes from k ^ modes. In the case of 
three dimensional systems, macroscopic occupancy of a 
single particle state is absent in one dimension, so k 7^ 
momentum excitations have much more significant effect 
in one dimensional dynamics. 

2. k 7^ contribution 

The exact evaluation of spin moments (^(Sf) m ^ for 

k is possible through the following trick. Con- 
sider the annihilation operator "f s . k (t) for the state 
such that 7s,fe(i) IV'C*)) = 0- If we write the 
operators </> s (r) in terms of J s ,k(t) and Ts fe (^) ' t nen 
k =/= part of the mth moment schematically takes 

the form ((Sf)™) ~ (exp(i£^ C., k 1,j* + C^J,*)) 
(Here and in the following, we drop the time de- 
pendence of 7 s ,fc(t) from the notation). Using the 
property e^-» \i/>{t)) = (1 + J s , k + • • • ) |V>(*)) = 
| -0 C*) ) an d tne identity e j4+s = e A e B e~^ A ' B ^ where 
[A, B] is a c-number, we can evaluate mth moments 
as ((Sf) m ) ~ (^ iC t^U e -hT, k ^o\CsM 2 e iC., hl s,^ = 

It is straightforward to check that 7 Sj fc operator is given 
by the linear combination of b Sj k and bl _ k as follows, 



is A*) 



-2W k e l 



-2ff t e'" 



(15) 

7^/c and 7] fc obey a canonical commutation relation 
[7s,fc) 7I J = !• In terms of these j s ,k, the expression 
of 4> Si k(t) becomes 



operator of our state at time t, we immediately conclude 
that (\^ k {t)\ 2 ) /L= \C Sik (t)\ 2 - 

Using the technique described above, mth moment of 
Sf becomes (we include both k = and k 7^ contribu- 
tions in the expression below) 



1/2 



W)\ [ I S»{r)dr m)) 
-1/2 J 



n 




\dn y: 



W.(r ( )-e) 



' 1 - .,=±1 

1/2 



' ' t f*l j£ 7 

^ j-j / pdTi ^ e (i( Sl 5 ( ri )+...+ Sm 4(r m )) 
{ Si =±l}i=l"^~ Z / 2 2 



e n 

{ Si =±l} i=l 



{s<,n}/t{si,n}\* 



(17) 



where £*' p '> - ^ ^f )?) { Sl e^ +... + s m e ikr ™). s t 
takes either the value 1 or —1, and E{ s } sums over all 
possible set of values. Note that L is the total system 
size and I is the integration range. 



3. Full Distribution Function 



Calculation of the full distribution functions from mo- 



ments in Eq. ( 17 ) is studied by the techniques introduced 



in Ref. [18j through mapping to the statistics of random 
surfaces. In this subsection, we provide the details of the 
calculation. 

Eq. (171 is simplified if the integrations for each 



can be independently carried out. This is not possible in 
Eq. (17) because e lkri and e lkrj for i 7^ j are coupled in 



{si,n} 



(l<<r ! ) + (*<fc ,r,} ) • To unentan- 
gle this, we introduce Hubbard-Stratonovich transforma- 



tion, e 2 



1 

/2rr 



e 2 e lxX , for example, 



1 

7l 



a 



s.k 



- e ic 3 \k\t _ 2W k e~ ic ^ t 
2\k\K s L ^l-4\W k \ 2 



.(16) 



C Sik (t) measures the fluctuation, or variance, of phase 
in the fcth mode at time t, given by ^|0 s ,fe(t)| 2 ^> = 

\<Ps.k{t)4 l s-k{t) )■ Indeed, since j S}k is the annihilation 



j(H«(fl'f' } )) 



d\is 



s /2 



We apply a similar transformation for Im^ SJ t. This re- 



moves the cross term between e %kTi and 



Akr; 



for i ^ j 



and allows us to independently integrate over r^'s. Asso- 
ciated with each transformation, we introduce auxiliary 
variables X lsk for Re{£ Syk ), A 2sfe for Im(£ Sife ). Then, mth 
moment becomes 
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m)\ [ 

En 



1-1/2 



-1/2 



S 9 (r)drj \m) = 

-(A?. h +Ag. t )/2 ^1»fc dX2sk 
V2n V2n 



n 



1/2 



-i/2 



exp 




r 



lsfc 



Re(e,) + A 2sfc Im(^V) - , 



where we introduced £j r ^ = < ^ a .^ t L T e »fcn_ Sum- 
mation over {s, = ±1} can now be carried out. Fur- 
thermore, we introduce a new variables X rs k and Xe s k, 
and replace Ai s fc and X2 S k through the relation X rs k — 
V X isk+ X Lk and cos(A esfc ) = X 2sk / V / Xj sk + X 2 sk . 
These operations result in the simplified expression, 



m)\ 

X Xrsk& 

where 



(1/2 \ m 

L/ {r)dr ) m) = ni/ 



rfA 



i/2 > 

/> / dr cos [x(r, {Ajak}) -0] 

'-i/2 , 



ask 



(18) 



x(r,{A isfc })=^ 



I'M 2 ) 



A rsfc sin(fcr + A esfe ), (19) 



'Js,k\ 



it sin 2 (c s \k\t) + al cos 2 (c s \k\t) 



2\k\K s 



Oik 



c s irt 

~k7 



VP 



(fc = 0), 



(MO), 



(20) 



with a k = The integration over X rsk and X Ss k in 



Eq.([18| extends from to oo and from — tt to 7r, respec- 
tively. 



Comparing the expression in Eq.(18l and the implicit 
definition of a distribution function in Eq.Q, it is easy 
to identify the distribution function as 



X rsk e~ x ^ 2 dX 



rsk 



x6\a-pj 1 ^ drcos{ X {r,{X jsk })-6]J . (21) 

This function can be numerically evaluated through 
Monte Carlo method with weight A,. s fee _A ™ fc / 2 for X rsk 
and equal unity weight for Xg sk . 

While we have assumed that the chemical potentials of 
spin-up and spin-down atoms are the same in the absence 
of interactions by going to the rotating frame, it is easy 
to obtain the expression for distribution functions in the 



lab frame. The energy difference E between spin-up and 
spin-down atoms results in the rotation of the spin in the 
x — y plane at a constant angular velocity E. Therefore, 
the distribution in the lab frame is obtained by replacing 
6 ->9 + Et in Eq. ((21). 

In this section, we have focused on the distribution 
function of spins in x — y plane, but it is also possible 
to obtain the distribution function of z component of the 
spin, and we present the result in the Appendix |A2| 



4- Joint Distribution Function 

From the expression for the spin operators in Eq. (|9|, 
we observe that the spin operators for the x and y di- 
rections commute in the low energy description. This 
is because spin operators in Tomonaga-Luttinger the- 
ory are coarse-grained over ~ p£ s particles, and since 
for weak interactions p£ s ^> 1, the uncertainty of mea- 
surements coming from non-commutativity of Sf and 
Sf becomes suppressed. The possibility of simultaneous 
measurements of spin x and y operators implies the ex- 
istence of joint distribution functions P; X,y (a,/3), where 
P^' v (a, f3)dad(3 is the probability that the simultaneous 
measurements of Sf and Sf give the values between a 
and a + da, and (3 and f3 + d(3, respectively. Here we pro- 
vide the expression for Pj c ' y (a, j3) and proves that this is 
indeed the unique solution. 

The joint distribution function Pf (a, /3) is given by 
the following expression 



n 



dX 



O.i k 



2 TV 



XrskS 



>x /2 dX r 



x 6 [a + i/3 — p 



1/2 



dre ix(r,{^s k }) 



1/2 



(22) 



where the expression for x(r, {Xj S k}) is given in Eq.( 19 1. 
To prove it, we first show that E q.(p2| ) reproduces the 
distribution function Pf(a) in Eq.(21| for all 8. Then, 



we show that a function with this property is unique, and 



therefore the expression in Eq.(21 1 is necessarily the joint 
distribution function. 

Given a joint distribution function Pf' y (a, f3), we can 
determine the distribution function Pf (7) of a spin point- 
ing in the direction (cos#, sin#, 0). Consider the spin S in 
the x—y plane with S — (a, ft, 0) whose probability distri- 
bution is given by Pj°' v (a, /?). The projection of the spin 
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S onto the axis pointing in the direction (cos 9, sin 9, 0) 
is given by |5| cos(0 - 9) where |5| = yfa 2 + (3 2 is the 
magnitude of spin and 4> is the angle Arg(a + i(3). After 
a simple algebra, we find \S\ cos((f> — 9) — acos9 + f3sin9. 
Then given a spin S = (a, /?, 0), if one measures the spin 
along the direction (cos 0, sin 0,0), the measurement re- 
sult gives 7 if and only if 7 = a cos + /3 sin 9. From 
this consideration, the probability distribution that the 
measurement along the direction (cos#, sin9, 0) gives the 
value 7 is given by 



J?( 7 ) = J dadpP?' v 



(a,/3)£(7- 



a cos ( 



/3 sin 0). (23) 



Now, if we plug in the expression of Eg. ( 22 ) in Eq. 
we see that Pf{i) agrees with Eq. (21 1 for all 9. 



(231, 



Now we prove the uniqueness of a function with the 
above property, i.e. a function that reproduces Eq. (21) 
through the relation Eq. ( 23 ) . Suppose you have another 



distribution P?' y (a, (3) that satisfies Eq.(23) for all 9. We 



define Q(a,/3) = Pf' v {a,0) - Pf' v {a,0). Our goal is to 
show that Q(a, (3) must be equal to zero. By definition, 
we have the equality 







J dad/3Q(a, 0) 8 (7 — a cos 6 — /3 sin 



6), (24) 



for all 9 and 7. If we take the Fourier transform of both 
sides of Eq. (24 1 in terms of 7, we obtain 



d<y / dad/3Q(a, /3)8(j - a cos 9 - (3 sin6»)e™ 7 



= / ddQ(a)e v 







In the last line, we defined w = w(cos9 1 sm9) and a = 
(a, (3). Notice that this equation holds for any w. Then 
this last expression is just like (two-dimensional) Fourier 
transform of Q. By taking the inverse fourier transform 
of the last expression in terms of w, we find 



= 



dw / daQ(d)e 



iw- (a — a f ) 



Q(a>), 



thereby proving the uniqueness of the joint distribution 

Pt y (a,{3). 




FIG. 6: The illustration of the dynamics for each harmonic 
oscillator mode, described by the Hamiltonian Eq.(|8|. The 
initial state contains a large fluctuation of density difference 
n s ,k given by (n Sj fen Si _fe) = ryp/2(see Eq.(lOl), and its con- 
jugate variable, the phase difference <p s ,k> has a small fluc- 
tuations. In the subsequent dynamics, such squeezed state 
evolves and energy oscillates between the fluctuations of the 
density difference and phase difference. 



Interpretation of the distribution dynamics 



The form of the distribution function in Eq. ( 22 ) 



encapsulates the interpretation in terms of dynamics 
originating from spin waves explained in Sec [TTJ. Here 
e ijc(i",{% S fc}) represents the spin direction at coordinate r, 
where the x — y plane of the spin component is taken 
to be a complex plane. Then Eq.(22) suggests that for 
a given instance of the set {Aj- S fc}, (Sf + iSf) is sim- 
ply the sum of the local spin directions e^^'^i** 1 ^ over 
the integration length I. The local spin direction at po- 
sition r are determined by the phase %(r, {tj s k}), which 
receives contributions from each spin wave of momentum 



From the joint distribution function in Eq.( 22 ), one can k with strength A k (t) = A 



also obtain other distributions, such as the distribution 

□1 
I 



(t) 1 2 > . Spin waves 



P^ij) of the square of the transverse spin magnitude, 



(Sy-(i)) , which is given by 

/CO 
dadf3P^(a,f3)S( 1 -a 2 -(3 2 ) 
-OO 

k j ° 

2\ 



xS 7 



8sk / \ 
27r Jo 

1/2 

-1/2 



(25) 



with momenta k rotate the spins as sin(fcr + Xg s k) (see 
the expression of x( r A^jsk}) in Eq.(19)). The rotation 
strength Ak(t) oc \ rs k has the distribution \ rs k^ > ' rsk ^ 2 , 
which represents the quantum fluctuation of the spins. 
On the other hand, Xg s k is distributed uniformly between 
— 7r and 7T. 

The dynamics of phase fluctuations ^|^ s ,fc(i)| 2 ^ can 

in fact be easily understood by considering the Hamilto- 
nian given by Eq.([8| as a harmonic oscillator for each 
fc(Fig. [6]). We first note that the initial state has a 
large fluctuation of density h s ^ because the initial 7r/2 
pulse prepares the state in the (almost) eigenstate of 
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K s = 30 K s = 25 K s = 20 



time=0 time =7 time =14 time =0 time =7 time =14 time =0 time =7 time =14 




-10 1-10 1-10 1-10 1-10 1-10 1-10 1-10 1-10 1 



FIG. 7: The dynamics of the joint distributions for L/£ s — 200, £ s = 40 and various spin Luttinger parameters K s — 30, 25, and 
20. Here axes are scaled such that the maximum value of a and j3 are 1. Smaller value of K s enhances the spin fluctuations, 
leading to a stronger diffusion and decay. Time is measured in units of £ s /c s . 



S x = pcos((j) s .k) with a small fluctuation of <^ s ,fe, the con- 
jugate variable of h s .k- The fluctuation of h s ± is given by 
(^s,fc"s,-fc) = 7 7P/2 (see Eq.(lOl). Because of this large 
fluctuation in the density, almost all the energy of the 
initial state is stored in the interaction term |n S! fc| 2 in 
Eq. Therefore, the total energy of each harmonic os- 
cillator can be estimated as 7 ^jf 1 - During the dynamics 
dictated by the harmonic oscillator Hamiltonian, this en- 
ergy oscillates between the density fluctuations and phase 
fluctuations in a sinusoidal fashion, see Fig. [6| In the 
dephased limit of the dynamics, approximately equal en- 
ergy of the system is distributed to the phase and den- 
sity fluctuations, and from the conservation of energy, 
we conclude that the characteristic magnitude of phase 
fluctuation is given by {\4>s,k(t)\ 2 ) ~ ak^u 2 ■ ^ ucn V^ 2 
dependence of ^|</>s,fcft)| 2 ^) agrees with the more rigorous 
result in Eq.(20). Therefore the spin fluctuations domi- 



nantly come from spin waves with long wavelengths, as 
we have stated in Section [TTJ Moreover, the weak depen- 
dence of spin dynamics on high momenta contributions 
justifies the use of Tomonaga-Luttinger theory for de- 
scribing the dynamics. We will more carefully analyze 
the dependence of distributions on the high momentum 
cutoff in Sec lIVDl 

From the simple argument above, it is also clear that 
the spin fluctuations coming from spin waves with mo- 
menta k have the time scales associated with the har- 
monic oscillators given by jj^r- Again, this rough ar- 
gument agrees with the more rigorous result presented 
in Eq.(20). Therefore, the fast dynamics is dominated 



by spin waves with high momenta and slow dynamics is 
dominated by low momenta. These considerations lead 
to the illustrative picture of Fig(2] Furthermore, this im- 
plies that the dynamics of the magnitude of spin (j3i~ 

reaches a steady state around the time since spin 
waves with wavelength longer than / do not affect the 



magnitude. This should be contrasted with the evolu- 
tion of the x component of spin which, in principle, keeps 
evolving until the time scale of ~ ^- (see Figji I . 

The strength of interactions and correlations are as- 
sociated with Luttinger parameter, K s . K s influences 
the spin fluctuations (\4> Sy k(t)\ 2 ) at all wavelength, and 
(|</ ) s,fe(i)| 2 ) depends on K s as \jK 2 s for a fixed density. 
As is expected, in the limit of the weak interaction corre- 
sponding to large K s , the amplitude of spin fluctuation 
decreases. In Fig. [7| we have plotted the time evolution 
of the joint distributions for L/£ s = 200, l/£ s = 40, and 
K s = 20, 25 and 30. For larger K s , we see that the spin 
fluctuations get quickly suppressed. 



C. Dynamics of the expectation value of the 
magnitude of spin /^§i~(t)^ 

In order to illustrate the dynamics of the Ramsey se- 
quence further, it is helpful to study the dynamics of the 
expectation value of the squared transverse magnitude, 



by((^(t))' 



given 

In Fig 8 we plot the evolution of 
K s = 20, L/£ s = 200 and l/£ 
ted (Sf(t)) along with 



with 



20, 30, 40. We also plot- 

2\ 



with the same 



parameters. It is easy to verify that (S[(t)) is indepen- 
dent of integration length As we have discussed 

in the previous section, J^Jjtf-^b^ ^ reach the steady 

states at the time scale of -J— with finite values, while 

4c s ' 

Sf(t) ) keeps decaying for much longer time. 
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FIG. 8: The dynamics of the average value of the magnitude 
of spins, a / ^S^-(t)^ V and the average of the x component 



of spins ^Sf(t)^. Here y axis is scaled such that the initial 

values take the maximum value of 1. Here we took L/£ s = 
200, K s = 20 and the integration lengths l/£_ s = 20,30,40. 
The magnitude of spins decays only due to the spin waves 
with wavelengths shorter than the integration length I, and 
the decay of the magnitude stops around the time scale of 
~ On the other hand, all spin waves contribute to the 

evolution of the of the x component of magnetization, which 
keeps decaying [58]. 



Si~(t)j j for sufficiently large integration length I 

attains the value which corresponds to the one expected 
from the randomization of spin patches of size £ s . At low 
energies, spins within the length ~ £ s are aligned in the 
same direction, but spin waves can randomize the direc- 
tion of the spin for each of Z/£ s patches. Since the magni- 
tude of spin within £ s is if the patches are completely 
randomized, the result of the random walk predicts that 

' (£,sP) 2 (l/£,s)- We will see below that, due 
to the properties of correlations in one dimension, the in- 
tegrated magnitude of spin y ^ (^) ^ never attains 

this form, albeit a similar expression is obtained (see 
Eq.(26l). Moreover, we identify the integration length 




I which separates the "spin diffusion regime" and the 
" spin decay regime" by finding the decaying length scale 



for 



, .(*) 

The results for the long time limit of 



can 



be analytically computed. Following similar steps leading 
to Eq.dltb, we find 



drpe 



Here &l } = \C., k \(e ikn 



a ikr 2 



). We introduce dimen- 



sionless variables r[ = ri/l,k' = kl and the integration 
over k in the exponent can be carried out as 



/ 

2 

7T 



2nl 



1 



dk' [ - ! -cos 2 (|fe|c s t) 



ir 2 pr]l 
k' 2 K 2 



sin 2 (|fc|c s i) 



x sin 



I A;' I 



kl 



2 P Trlri 2K S 1 



dy 



sin 2 (y) 



y- 



It is interesting to ask if the long time limit of Since J Q dy 



In the second line, we approximated cos 2 (|fc|c s i) rs 
sin 2 (|£;|c s i) w 1/2, which is appropriate for long time. 
In the last line, we extended the upper limit of the in- 
tegration for the second term to oo and the lower limit 
to 0. The former is justified because we know that high 
momentum contribution is suppressed by 1/k 2 , and the 
latter is justfied because we also know low momenta ex- 
citations with wavelengths larger than I do not affect Sj . 

7r/2, we find, in the long time limit, 



00 j,, gig fa) _ 



(5^(t = oo)) 2 \ / 0)} 



(1 - exp(-Zo)) 



(26) 



where we expressed the result as a ratio of the asymptotic 
value and the value at shortest time scale of the theory 
given by t ~ 1/ p. Iq = IS ^ ne dimensionless integra- 
tion length that controls the value of 
long time limit. As soon as Iq becomes larger than 1, the 



long time value of ( ( 



intcgration limit, oc 2 



quickly approach the long 

2 ' 



Therefore, Iq ~ 1 



separates the " spin diffusion regime" and the " spin decay 
regime." 

An intuition behind the expression for Iq can be ex- 
plained through the following heuristic argument. The 
system enters the spin decay regime when the spins 
within the integration length I rotates by 2-7T across I. 
The angle difference between the spins at r — and 
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I in the long time limit is roughly given by A% 



7Z Sfc ^rsk \j( \<f>s,k\ 2 ) sin(fci) where (\<j> 3> k\ 



is the characteristic magnitude of \\4> s ,k{t)\ 2 / m Eq.(pO|, 
which is given by the half of the maximum magnitude of 
(\4> s ,k(~t)\ 2 ^' Now the expectation of magnitude ((Ax) 2 ) 
over the quantum fluctuations represented by X rs k can 
be computed, and it yields ((Ax) 2 ) « §§^. When 

\J ((Ax) 2 ) becomes of the order of 1, the system enters 
the spin decay regime. This estimate gives the boundary 
between the two regimes Iq = w 1 apart from an 

unimportant numerical factor. 



It is notable that the Eq. ( 26 ) approaches the random 
walk behavior oc (£ s p) (£/£«) very slowly, i.e. in an al- 
gebraic fashion. Therefore, even in the steady state, the 
system retains a strong correlation among spins. More- 
over, Eq. (26) in the limit of Iq — > oo is not just the 



random walk value, but is proportional to K s , which 
measures the strength of fluctuations. 

The calculation above shows that the spin diffusion 
regime and the spin decay regime are separated at the 
integration length scale of I « This length scale 

is nothing but the correlation length of spins in the long 
time limit. The calculation of the spin correlation length, 
for example, between S x (r) and S x (r') can be done sim- 
ilarly to the calculation of /^jSj 1 "^ ). The resull in tlic 
long time limit is 



S x (r)S x (r'] 



-r'\/t 



(27) 



where C 



-k c / (47rp7/) 



is a small reduction of the spins 



due to the contributions from high energy sector. Thus, 
one expects qualitatively different behaviors of distribu- 
tion functions for integration lengths I < I and I > I. 



that (here we ignore k = contribution) 



S?) = 



1/2 



1/2 



=pl exp 



^dx (e ilt>( - x) -r-e-^*^ 



J2\c s , k \ 2 = [ c 



k " /cos 2 (|fc|c s i) npri . ,,,,, 
dk [ } l 1 ' + , , "J rn sin 2 (|fc|c s i) 



Airpri 4k 2 K 2 

k c pc s ti] 



2^pr 1 {2K s /nf 



(28) 



where in the last line, we took the long time limit 
t £ s /c s [68]. In this limit, only the first term in 
Eq.(28) depends on the cutoff k c , and moreover, the 
cutoff dependence is independent of time. The effect 

is to reduce the value of (Sf^ through the multiplica- 
tion of a number close to one in the weakly interacting 
limit. For example, if we take k c = 27r/£ c , then the 
cut-off dependent term reduces the value by multiplying 
eK P (- 3 ^)«e-V(«.)«l. 

In a similar fashion, higher moments of spin opera- 
tors can be shown to have a weak dependence on the 
cutoff momentum k c , as long as the integration length 
is much larger than the healing length, l/£ s ^> 1. In 
this limit, m moments of, for example, Sf is reduced by 
-m-^- I. Therefore, the full distribution function 



exp 



V 4 ^pv J 



is simply reduced by the multiplication of a number close 
to one exp ^— jf^) ~ e^ 1 ^ 4K ^ w 1 in the weakly in- 
teracting regime. This gives the self-consistency check of 
our results in Sec IIVBI 



V. DYNAMICS OF FULL DISTRIBUTION 
FUNCTION IN THE PRESENCE OF MIXING 
BETWEEN SPIN AND CHARGE DEGREES OF 
FREEDOM 



D. Momentum cut-off dependence 

The description of dynamics presented above uses the 
low energy effective theory. In order to confirm the self- 
consistency of our approach, we check that the distri- 
butions of spins are not strongly affected by high energy 
physics, i.e. they weakly depend on high momentum cut- 
off. We have seen an indication that this is indeed the 
case through the weak fluctuat ions of phases for large k, 

\4>s,k\) oc 1/fc 2 , in Sec 



IVB5 



First of all, we analyze the high momentum cut-off 
k c ~ 2-7r/£ s dependence of the average value of Sf. From 



In this section, we extend the analysis in Sec |IV| to 
a more general case, in which spin and charge degrees 
of freedom mix. We will see that the distribution func- 
tions even for this more general case have essentially the 
same structure as in Eq. (21 ), and are described by spin 



the discussion in Sec IV B[ it is straightforward to obtain 



waves with fluctuations whose amplitude is determined 
by the fluctuations of phase ^l^s.fcl 2 ^- One important 
difference from the decoupling case is the dependence of 
spin distributions on the initial temperature of the sys- 
tem. The thermal excitations are present in the charge 
degrees of freedom in the initial state, and such thermal 

fluctuations increase the value of (^\4> s ,k\ 2 ^j through the 
coupling between spin and charge during the evolution. 
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Hamiltonian and initial state 



In a generic system of two component bosons in one 
dimension, spin and charge degrees of freedom couple 
through the mixing Hamiltonian in Eq.Q. Yet, Hamil- 
tonian in Eq.Q is still quadratic and it can be diagonal- 
ized. We define new operators <f>i, <f>%, hi, h% by 



n 

4>2 

h-2 



COS K sin K 
— sin k cos k 



COS K sin K 
— sin k cos k 



(29) 
(30) 



Mixing angle k and scaling parameter s c are chosen so 
that the Hamiltonian is written in the following diagonal 
form 



H = 
Hi = 

Ho = 



H x +U 2 

L/2 

dr 

-L/2 2m i 
L/2 



(31) 



P 



dr- — 
L/2 2m 2 



(V^ 2 ) 2 +92(h 2 f 



Explicitly, n and s c are given by 



2m c 



- + 9s9 



mix 



2m a 



K 



s c 9c ~ 9s 



2m c 



tan k = 



,c 2m, 



-«o ± \Ao + 4 



fJn 



where ± in the expression of tan k is + when k > 
and — when kq < 0. We defined k such that k — 
corresponds to decoupling of charge and spin, i.e. to 
9mix = and g^ mix = in Eq. (m| . Parameters gi, g 2 , §^ 
and t£— are given by 



2m 2 



9l 

<)2 



2mx 

P 
2m 2 



s c9c + ^/~^c tan t^Qmix i 

9 s ~ vSctanK^ m j X , 
P 1 tanK- mix 



2m r s r 



2m, 



tan re- 



in the weakly interacting systems which we study in this 
paper, Luttinger parameters Ki and sound velocities Cj 
are determined for each Hamiltonian Hi, i = J [, 4,, c, s, 1, 2 
through 



K,=n 



P 



2m i g l 




(32) 



At finite temperature, the state before the first ir/2 
pulse contains excitations, and these excitations are car- 
ried over to the charge degrees of freedom after the pulse. 



Pulse only acts on the spin degrees of freedom, and the 
local sum density of spin-up and down is left untouched 
as long as the pulse is applied in a short time compared to 
the inverse of typical excitation energies, (3 = l/(ksT). 
In other words, the local density fluctuation of spin- 
up, h^(r), before 7r/2 pulse is converted to the sum of 
the local density fluctuation of spin-up and spin-down, 
tv\{t) 4- h^(r) after 7r/2 pulse. In this strong pulse limit, 
then, the distribution of hf(r) before n/2 pulse is the 
same as the distribution of n^(r) + n±(r) after tt/2 pulse. 

The distribution of the local density for spin-up atoms 
before 7r/2 pulse is determined by the density matrix for 
spin-up given by e~^ H ^, where in the weak interaction 
regime we have (see Eq.Q) 

-L/2 



H^ = 



dr 



L/2 



2p_ 

2m^ 



(VM r )) +9n(M r )Y 



Then, the density matrix which produces the distribution 
of hf(r) -\-h±(r) required above is given by e~$ H ^ where 



H \o = !^L/2 dr 



L/2 



2p 
2 m 



-{(V0 t (r) + V^(r))/2} 2 , 



+g n (h t (r) +h i (r)) 2 



f -f-iUr 



L/2 ' 



^(V0 C ) 2 



(33) 



where K<^ = 



and c c | 



m^ 



The initial state for spins is determined by the n/2 
pulse, and we obtained the state in Eq.(12). Then, the 



complete initial density matrix after the first n/2 pulse 
is given by 

Po = l^o) (H ® e-^t/Tr . (34) 

This density matrix evolves in time as pit) — 
e~ ltH poe ltH . Since we assume that the preparation of 
the initial state is done through a strong, short pulse, 
the spin and charge degrees of freedom are unentangled 
in the initial state. 



B. Time evolutions of operators 

In order to calculate the distribution function of 5?, 
we again start from the calculation of mth moments, 

Tr (j)(t) {sf^j ^ . Evaluation of moments can be done 

through a similar technique used in Sec |IV Bj 

In the following, we describe convenient, time- 
dependent operators 7 s ,fc(£) and J c ,k(t) used to evaluate 

spin operators such as e^ s fc . The first operator resides in 
the spin sector and it is again the annihilation operator 
of the initial spin state such that Tr7 s> / c (0)/5o = 0. This 
operator is given in Eq. (15), which is 

~ itH -fsAoy tH , 

/ 1 -2W k 

7l- fe (0) 
7.,k(0) 



ls,k{t) = 



yJl-4\W k \ 2 y/l-4\W k \ 2 

-2W k 1 
y/l~4\W k \ 2 y/l-4\W k \ 2 



,-k 



Js,k 



17 



with 2W h 



^ and a k 



\ k \ K * as before. The second 

7rpr] 



operator is the operator of charge degrees of freedom, and 
it is given by 

7c,*(f) = e- itH lcjk (0)e itH , 
7c,fc(0) = b c -f :k - 

where 6 c ^.fc is an annihilation operator for the elementary 
excitations in H c ^. Since 7 s ,fe(t) and 7 c ,fe(t) commute at 
t = 0, they commute at any time t. We will drop the 
time dependence of j a ,k(t) m the notation from now on. 
From the expression of initial density matrix po in Eq. 



( 34 1 , it is easy to check that the density matrix at time t 
given by p(t) = e~ ltH poe ttH can be written as the tensor 



product of the density matrix of operators j s ,k(i) an d 
that of 7 c ,fc(t). This is because po is a tensor product of 
the density matrices of 7 S) fc(i = 0) and that of J c ,k (t = 0). 
This structure of the density matrices at time t allows the 
independent evaluation of "f s ,k(t) and "f c ,k{t) operators, 
and it is advantageous to express spin operators in terms 
of these operators. 

As we show in the Appendix B we can write <fi S) k in 
terms of 7 c ,fe(i) and 7 S) fe(i) as follows. 



'L 



*>s,k = C* fc 7 s f _ k + C S:k j s , k + ^*,)fe7c + C c ,fe7c,fc, 

(35) 

where explicit expression of C s ^ and C c ^ are given by 



Cs,fc 



2Lpr) 



{cos 2 9 cos(c2|fc|£) + sin 2 0cos(ci|fc|i)} — i 



K s ( cos 2 #sin(c2|fc|£) sin 2 0sin(ci|/c|i) 



Q'A: 



C c k = cos 9 sin ( 



2L\k\s r K, 



i {cos(ci|fc|i) — cos(c2\k\t)} — 



K 2 K x 
sin(c2|fc|t) sin(ci|fc|t) 



Ki 



(36) 



where K c ^ = K c ^ / ^fs 



of 7a, k with a = s, c as follows 



Using Eq.(35 ), we find an expression for (Sf ) m in terms 




where £, a ,k = (J2T s % eirik )C a .k- In the following, we sep- 
arately evaluate three contributions; k = component 
given by e l( ^ B i)4',o/vL- the charge component of fc 7^ 
given by e I ^ fc ^ ^ c fc7c ' c+ ^ c ' fc7c ' fc ' ) ; the spin component of 
k =/= given by e l ^fc^o(C,fc7], fc +« s ,fc7 s ,fc)_ 



details to the Appendix [Cj and the result is 




contribution 



The initial state of k — spin sector in Eq.(12| as 



well as that of the charge sector in Eq. (331 both have 
a Gaussian form so that calculation of the trace Tr 
je^E.sO^.o/VE^)! i s straightforward. We leave the 



2. k 7^ 0, spin sector 



This calculation is analogous to Eq.( 17) and the result 



can be directly read off from Eq.(17l, and it is 



eXP ~lj2^tkis,k ■ 

V Mo / 

(38) 
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3. k 7^ 0, charge sector 
We first rewrite the density matrix at time t as 

where Af c is normalization given by Af c = 
Tre -Pc c ,J2 k ^\k h l k (th c , k (t) = H k ^ -1/X k with 
X k = er/^tlfel - 1. 

Then the trace of I Sf J for fc ^ spin sector is 

^ e l E fc; io(?c,fc7 c t ifc +«c,fc7c,fc)^ _ ^g-j 
fe^O 



We can evaluate this by taking the trace in the 
basis of normalized coherent states \a k ) such that 
7c,fe l^fc) = a k \<*k)- The use of the identity 1 = 
^ J cPctk \otk) (ak\ as well as of an important equality 
e va 1 a _, e (e -i)a f a . ^^SJ , where : : is a normal ordering 
of 0, leads to 



e 4 E fe ^o(C,fcT c t ,fc+?c,fc7c,fc) 



i (ic,fcTc t ,fc+?<=.fc7c,fc)„-/3|fc|c ct 7e lfc 7c,fc 



l , , i r 

=JfI[- £ak < a *l e 

' k^O J 



MO 
1 i — r 1 

= /7 



MO 



1 1+ ,-g;*;°«t 



^. .FWZ distribution function 
We can summarize the results above as 



(Sf )") = E ft / ex P f -\H C.kUk I exp I -\ E Mcfc^^c* ] exp (~ ( 



-i(£ 4 *)0 



(40) 



Here, M c j, = ^ +e _ )3|fc|eef . As before, we introduce 

the auxiliary variables to separate spatial integrations 
over r 4 . We can combine £* fc £ s ,fc and M c ^* ck ^ c , k so 
that we only need to introduce three sets of variables, 
Ai, s ,fc, ^2,s,k, Ao for Hubbard-Stratonovich transforma- 
tion. Summing over {s^} simplifies the result, leading to 



the following expression for the full distribution function 

if («) - Ilk ^ JT W^-/ 2 dA rsfe , 
x c) (a - p cos [x(r, {A jsfe }) - 6>]) 

(i^, fc i 2 )/L= ic s , fc i 2 + ^ i :;i:; ic c , fc i 2 , 
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a) Qmix = 0.1g s b) 9mix = 0.3c7 s 

time=0 time =7 time =14 time =0 time =7 time =14 




FIG. 9: Time evolution of the joint distribution function P x ' y (a,/3) for the system size L/£ s = 400, the spin Luttinger 
parameter K s = 20 and integration length Z//£ s = 20 in the presence of mixing between the spin and charge modes. For a), 
the interaction strength ratio is taken to be g c : g a : g m ix = 1:1: 0.1, and for b), g c : g 3 : g m ix = 1:1: 0.3. Time is measured 
in units of £ s /c s where c s is the spin sound wave velocity. Here axes are scaled such that the maximum value of a and /? are 
1. With increasing strength of mixing, the large initial temperature affects the spin dynamics at earlier time more strongly. 



The last line can be confirmed by directly computing 
2\ + v,„ ;„ v„ 1^, The ex pres- 

As before, the 



\4>s,k\ ji using the expression in Eq. 

sion for (|^> s ,o| 2 ) is given by Eq. |37l 
joint distributions as well as the distributions of squared 
transverse magnitude can be obtained through the same 
procedure as in Sec |IV B 4} 



The spin distribution in the presence of mixing be- 
tween spin and charge degrees of freedom resembles 
the one in the absence of such mixing, and the only 
change is the additional contributions to phase fluctu- 
ations coming from the ther mal excitations, represented 
12 in Eq 



by \CcM 



(411. \C C 



(3GI 



is proportional 
Thus, for weak 



to sin 2 k as one can see from Eq. 
coupling of k ~ 0, the contribution is diminished by a 
factor of k 2 . 



In the experiment by Widera et ai, they used Rb 87 
in the presence of Feshbach resonance. They employed 
the theory which assumes the absence of mixing between 
spin and charge degrees of freedom to analyze the decay 
of the Ramsey fringes. The ratio of interaction strengths 
in their experiment can be roughly estimated as g c : g s : 
gmix ~ 3.66 : 0.34 : 0.06 which leads to the value of 
k « 2 x 10~ 2 . Therefore, the thermal contributions are 
diminished by about four order of magnitude and thus, 
their assumption of decoupling between spin and charge 
is justified. 



In Figj9j we have plotted the evolution of the joint 
distribution functions for different strength of the cou- 
pling g m i X at a relatively large initial temperature /c^T = 
0.4 x 2irc c ^/t; s where 2ttc c ^/£ s is approximately the high 
energy cut-off of Tomonaga-Luttinger theory. Here we 
took the system size L/£ s = 400, the Luttinger parame- 
ter K s = 20, integration length I = 20£ s . For Figj9] a), 
the ratio of interaction is taken to be g c : g s : g m ix = 1 : 
1 : 0.1, and for For Fig[9]b), g c : g s : g mlx = 1:1: 0.3. 
One can see that with increasing strength of mixing, the 
large initial temperature affects the spin dynamics at ear- 
lier time more strongly. For comparison, also see Fig. [4] 



VI. INTERFERENCE OF TWO 
ONE-DIMENSIONAL CONDENSATES 



A. Dynamics of interference pattern 



As we have described in Sec II B the full distribution of 



interference patterns can be studied in exactly the same 
way as we have studied the full distribution of spins in 
previous sections. In the following, we more formally 
describe the dynamics of split condensates. 

The low energy effective Hamiltonian of two quasi- 
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condensates after splitting is given by 



H = 
H L = 

H R = 



H L + H R 

L/2 

dr 

-L/2 
L/2 

dr 

-L/2 



(42) 



(V^(r)) 2 + §(M»0)' 



2m 2 



where we assumed weakly interacting bosons with a pos- 
sible density difference pr — pl between the two 
condensates. Here and in the following, we consider the 
rotating frame and ignore the chemical potential differ- 
ence g/2(p 2 R — p L ) between left and right condensates 
arising from interactions. 

The interference pattern measures the phase difference 
4>l — <j)R- We describe the system in terms of the "spin" 
variables that are the difference of left and right con- 
densates and " charge" variables that are the sum of the 
two. Using the variables cj) s = 4>r — <Pl, 4>c — 4>r + 4>l, 
n s = (n R — n/,)/2, fi c = (fiR + h^) /2 , we find the rlamil- 
tonian of the system to be 



H 

Hint 



H s + Hc + H int 

PR. + PL 



dx 
dx 
dx 



8to 

PR + PL 

8m 

PR - PL 
Am 



{d x (j>s? +gn 2 s 
(d x 4> c ) 2 +gn 2 c 
d x <t> c dx<i>. 



(43) 
(44) 

(45) 

(46) 



Therefore, when the splitting makes two identical quasi- 
condcnsates with equal density, "spin" and "charge" de- 
grees of freedom decouple and we can use a simpler the- 



ory derived in Sec IV B On the other hand, when the 
splitting makes two condensates with unequal densities, 
more general theory of Sec [V] needs to be employed. In 
any case, the full time evolution of the distributions of 
interference patterns can be obtained, which in principle 
can be compared with experiments. 

It is notable that the mixing of the "spin" and "charge" 
degrees of freedom for small density difference pr — pl is 
not "small," in the sense that the mixing angle n defined 
in section [V] takes the maximum value n/4. The spin 
decoupling in the limit of pr — pl — > is recovered not 
by taking k — > 0, but rather, by taking the time at which 
the effect of the coupling takes place to infinity. This 
is most explicitly shown in Eq. ( 36 ) where the charge 
contributions of fluctuations go to zero as c\ — > C2 which 
is attained in the limit pr — pl = 0. 



functions of the interference patterns between two inde- 
pendent, thermal quasi-condensates. This problem has 
been previously analyzed in theory [5U [7D] as well as in 
experiments [TH HT] . 

We consider the preparation of two independent one 
dimensional quasi-condensates. If they are prepared by 
cooling two independent quasi-condensates, the temper- 
ature of the left quasi-condensate Tj, and that of right 
quasi-condensate Tr are generically different. The den- 
sity matrix of the initial state is described by po — 
e -(hH L +0 R H R ) where p a = \/(k B T a ) with a = L,R. 

It is important to note that the constant shift of phase 
4>a 4>a + 9ac docs not change the energy of the sys- 
tem, so that for the average over thermal ensemble one 
has to integrate over 9 ac . Physically, this simply means 
that the phases of independent condensates are random. 
Then the only interesting distribution here is the distri- 
bution of the interference contrast [TH [18, 26, 69 given 
by 



1/2 



1/2 



-%4> s (r) 



dr 



(47) 



which corresponds to, in spin language, the squared 
transverse magnitude of the spin • The analysis 

of the evaluation of distributions in the density matrix 
of thermal equilibrium state in Sec [V] can be directly ex- 
tended to this case, and we obtain the distribution 



^(7)=n 



d\ 



2tt 



^rsk& 



x 6 



1/2 



-1/2 



(48) 



X(r, {Aj S fc}) 



\4>s,k\ 2 



\kk\ 2 ) 

L 



X rsk sin(/cr + Ag sfc ), 



I _|_ e -/3h\k\c L 



1_ e -h\k\c t 2\k\K L 
1 + e ~' 3 «l fe l c « 7r 



+ 



1 _ e -l3n.\k\c R 2\k\Kj- 



(49) 



where c a and K a ,a = L,R are the sound velocity and 
Luttinger parameters of left and right quasi-condensate. 



C. Prethermalization of interference patterns 



B. Interference patterns in equilibrium 

The techniques to calculate the full distribution func- 
tions presented in previous sections are directly applica- 
ble to also obtaining a simple form of the full distribution 



In Sec jll C] we gave a heuristic argument for prethre- 
malization phenomena, where the distribution of the in- 
terference contrast amplitudes of the two non- equilibrium 
quasi-condensates are given by that of two equilibrium 
quasi-condensates at some effective temperature T e g. We 
identified the effective temperature to be the energy 
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stored in each momentum mode. In Sec. IV B 5 we found 



this energy to be ^ff^, thus we conclude fcsTrff = 7r 4 c ^. p ' ? . 

In the following, we formally derive the result above, 
using the expressions of full distributions of interference 
patterns. The distribution of the inter fere nce contrast is 

2 \ given in Eq. ( 20 1 . In the long 



determined by (A (ft 

time limit, we can take siri^(c s |fc|i) ~ cos 2 (c s |/c|i) ~ 1/2. 
Moreover, since the interference contrast is most affected 
by the excitations with small wave vectors k with — 
< 1, we can approximate the expression as 



2\k\K s 2\k\K s 



(50) 



On the other hand, for two quasi-condensates in ther- 
mal equilibrium, the position of the interference peaks is 
again random. The inter fere nce contrast is determined 

by (\$8,k 



given in Eq.(49). Since the main contribu- 



tion to the fluctuation comes from low momenta, we ap- 



proximate e 



-P\k\c 



l — p\k\c. It is easy to check that the 



sound velocity and Luttinger parameters for each con- 
densate is related to those of the difference mode (see 
Eqs.( [44p6| )) as c L = c R = c s and K L = K R = 2K S . 
Thus we obtain 



\<t>. 



s,k\ 



P\k\c 8 2\k\K s 



(51) 



Now the crucial observation is that our closed form 
expressions for distributions of interference contrasts 
of both split quasi-condensates and thermal quasi- 
condensates are determined solely by ^|0 s ,fc| 2 ^), and they 

take precisely the same form in terms of ( |0 s ,/t| 2 ). More- 



over, the expressions given by Eq.(50l and Eq.(51) have 



the same dependence on wave vectors \k\. Therefore, 
the full distribution of interference contrast of split con- 
densates become indistinguishable from that of thermal 
condensates with temperature 



knT r 



AK„ ~ 2 ' 



(52) 



where the second equality holds for weakly interacting 
bosons and the chemical potential of one quasi-condesate 
is given by p, = gp. Thus, split one dimensional quasi- 
condensates indeed display the prethermalization phe- 
nomenon. 

In Fig . [ro| we plot the interference contrast P;" L (7) 
(see Eq.(|25[)) of split condensates in a steady state at 
time t — 60£ s /c s for Luttinger parameter K s — 20, sys- 
tem size L = 400£ s and two different integration lengths 
l/£ s — 30,40. Also we plot the interference contrast of 



the thermal quasi-condensates (see Eq.(48l) at tempera- 
ture ||*- for the same integration length. This tempera- 
ture corresponds to the effective temperature obtained in 
Eq. ( 52 1 . Indeed we see only a small difference between 
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FIG. 10: The distributions of interference contrast for steady 
states of split quasi-condensates and two thermal quasi- 
condensates. Here x axis is scaled such that the maximum 
value of interference contrast is 1. For the split condensates, 
we plot the distribution at time t = 60£ s /c a for Luttinger 
parameter K 3 — 20, system size L — 400£ s and two dif- 
ferent integration length l/£ s — 30,40. The thermal quasi- 
condensates are for temperature for the same integration 
length corresponding to the effective temperature obtained in 
Eq. pi. 



both integration lengths. The small difference comes 
from the approximations made in obtaining the expres- 
sions given by Eq.(50) and Eq.(51). 



In the previous paragraphs, we assumed that the split- 
ting prepares quasi-condensates with identical average 
densities. Here we briefly consider the case in which 
the splitting process prepares two quasi-condensates with 
slightly different densities. In this case, the temper- 
ature of the initial quasi-condensates affects the inter- 
ference contrast around the time scale of — 

h 



(ci— cr)tt 

whereas the prethermalized, long-time tran- 
where I is the integra- 



/PL-y/PR 

sient state is reached around - 4- 



the distributions of steady states and thermal states for 



tion length. 

These analytic arguments can be confirmed through 
numerical simulations. In Fig. |11[ we have plotted the 
evolution of the interference contrast C 2 for system size 
L = 500£ s , integration length I = 40£ s , and Luttinger pa- 
rameter K s — 20 with initial temperature corresponding 
to the chemical potential p. Here we consider a situ- 
ation where the density of left quasi-condensate is dif- 
ferent from that of right quasi-condensate by 20% such 
that pn = l.2p and p^ — 0.8p where pm^ is the aver- 
age density of the right (left) condensate, and 2p is the 
average density of the initial condensate before splitting. 
Also for comparison we have plotted the magnitude of 
interference contrast of two quasi-condensates in thermal 
states at temperature given by p. From the plot, one 
can observe the existence of quasi-steady state plateau 
after short time. Notice that the magnitude of C 2 in 
the steady state is larger than the value expected from 
thermalized states at the initial temperature. The subse- 
quent slow decrease of the interference contrast is due to 
the effect of temperature in the initial state coming from 
the small difference of the two quasi-condensates. Such 
development of the plateau at larger value of interference 
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Thermal condensates 



20 



40 



60 



Time 



FIG. 11: The evolution of the interference contrast C 2 for 
system size L = 500£ s , integration length I — 40£ s , and the ef- 
fective spin Luttinger parameter K s — 20 with initial temper- 
ature corresponding to the chemical potential p. Time is mea- 
sured in units of £ s /c s . Here we took the density pn = 1.2p 
and pL = 0.8/0. The magnitude of C 2 for two thermal quasi- 
condensates at temperature ksT = p is plotted as red dotted 
line for comparison. 



contrast than the one for equilibrium state of the initial 
temperature indicates the phenomenon of prethermaliza- 
tion. 



of interest to extend our analysis to higher spins [75], and 
analyze them, for example, in the presence of magnetic 
field[73]. Since there are more degrees of freedom in 
these systems, distributions might capture the tendency 
towards various phases such as ferromagnetic ordering. 
These questions will be analyzed in the future works. 
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APPENDICES 



Appendix A: Distribution function of the z 
component of spin 

In this appendix, we calculate the distribution function 
of Sf in the absence of the coupling between charge and 
spin. The extension to the case in which the charge and 
spin degrees of freedom mix is straightforward. 

It is convenient to evaluate the generating function 

^e^*J, instead of distribution function Pf{a). They 
are related by 



VII. CONCLUSION 

In this work, we have shown how noise captured by full 
distribution functions can be used to study the dynamics 
of many-body system in one dimension. The analytical 
results of joint distribution functions obtained in Sec |IV B] 
allow not only the simple understanding of distribution 
functions from spin-wave picture, but also an intuitive vi- 
sualization of the correlation in one dimensional system. 
Using this picture, we have also shown that the phenom- 
ena of prethermalization occur in one dimensional dy- 
namics. The thermal- like behaviors of the prethermalized 
state is revealed through the full distribution functions 
that contain information about the correlation functions 
of the arbitrary order. For the experimental demonstra- 
tion of such prethermalizations, see Gring et al[2A . 

The approach developed in this paper can be extended 
to other types of dynamics. While we focused on Ram- 
sey type dynamics or dynamics of interference patterns 
for a split quasi-condensate, we can also change different 
physical parameters to induce the dynamics. For exam- 
ple, it is straightforward to apply our study to the sudden 
change (quench) of interaction strength [4"Tjl I7T] . 

In this paper, we focused on distribution functions ob- 
tained from Tomonaga-Luttinger Hamiltonian Q. It is 



AS? 



e AQ P ; z (a)da. 



(Al) 



This equality can be checked by differentiating both sides 
by A and evaluating them at A = 0. This reproduces the 
implicit definition of Pf in Eq. ([2| . 

Analogous to the calculation of roth moment of Sf , we 
first express Sf in terms of j s k operators defined in Eq. 
PI 



Sf(r) 



J ~ 1 / 2 \k^0 

\k\K s e icMt + 2W k e- ic °W t 



1 



2ttL 



\/l-4|IU fe | : 



Then, we can apply the trick introduced in Section 
EH to obtain 



A 2 Jj>% dr ld r 2 ( £ fc , 1^,1 V^l-2UM 



e xs >) = e ' " 



exp < A 



E 



41 d„ 



AL ' ^ k 2 

fc#0 



sin 2 (7fc/2) 
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Then the following expression gives the distribution of 



( 



where 



: exp 



2WW IS, 



\ 



S, 



(A2) 



Sf 



pyi 2 



E 



4|4,fc| 
k 2 



sm 2 (lk/2). 



Appendix B: Expression for C a .k in the presence of 
mixing between charge and spin 

In this section, we derive the expression of C aik , a — s, c 
in Eq. (36). We first find the transformation from 
b s ,k,bcf,k to &i,fc,&2,fc- Then we relate bi jk ,b 2 ,k and 
7s,fc(i)j Jc,k{t)- Combining these two transformations, we 



obtain b St k in terms of ~/ s ,k(t)) 7c,k(t), leading to the ex- 
pression of (j) s , k in terms of 7 s ,fe(t), 7c,fc(£)- 
From the relations, 



-(6 



n it k 



2tt 



2\k\K, 



\k\K ifh1 



i4>i,k 



\k\Kj 
2ir 



riLk 



2\h\Ki' 



along with Eq. ( 29 ) , it is straightforward to obtain 



( b U,-* \ 

V b sM J 



= D 



( b U \ 

bi,k 
V b 2M ) 



where 




where K, 



Next, we relate "f a ,k(t), a = c, s operators to b±, b 2 . At 
t = 0, we have the relation between 7 a ,fc(0), a = c, s and 
b c ^.k and b S:k as described in SeclVl Since operators b c ^, k 
and b Syk are related to bi k and DjT fc through the matrix 
£> in Eq.( |Bl[ ), we can express 7 a ,&(0) as a linear combi- 
nations of bi] k and b 2}k . The time evolution of 7 a ,fc(0) is 
quite simple now because Hamiltonians are diagonal in 
the basis b\^ k and b 2k . These considerations lead to the 
relations 





-kit) \ 


Jc,k(t) 


7! 


-kit) 


V 7«,fc(*) / 


E k = F k D, 




(X 
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F k = 


V 




[ 00 v 



Ek 



( e -«i|fc|t & t _ k \ 

p -ic 2 \k\tjJ 
e ic 2 \k\t b 



<J 2 ,k 






-2H'i, 



Now define a matrix G(k) = DE k 1 so that 



b c ,k 



b\,- k 



G(fc) 

Then finally C^fc are given by 



lc,k 

ll,-k 

V ls,k ) 



a 



s . k 



Gr.k 



2L\k\K. 



(G 34 (k) - G 44 (fc)), 



2L\k\K. 



(G 32 (k) - G i2 {k)), (B2) 



where Gij(k) are the matrix elements of G{k). 

Appendix C: k — contribution in the presence of 
mixing between charge and spin 

In this appendix, we evaluate Trje'^i «*)**,o/>ff p (t)\ 
We first obtain the operator after time evo- 
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lution as 



iHt 



exp z 



Vl 



{Msfi + A'n sfi 



D 



B' 



' J c,0 



n c,0, 



(CI) 



Coefficients A, A' , B, B' can be found as follows, k = 
part of the Hamiltonian is given by H — ^§^n\ + 
^iijp (see Eq.jlTj). Using the commutation relation 

[ni,o,4>ifi] = -i, we have e im cj> ifi e- iHi = ljO + j^n ifi t. 
With the relation, </> Si o = sinK^o + cosk02,o, we obtain 

e «« e (£i -iHt = 



I sin K(0i i o + : ^j L ni,ot)+cos K(02,o + T;j«2,ot)} 



Now by transforming back to c, s basis through Eq.(29l, 
we find 



.4 
A' 
B 



1, 

sin K^——t 4- cos K-r—t, 
0, 



B = sin K cos k — — t. 

\Ki A 2 / 

Now that we know the operator after time-evolution Eq. 
(CI ), we evaluate it in the initial state. 

Initial state of the spin sector is \ip s ,k=o) m Eq. ( p~2| ) . 
Since this state is Gaussian, we have the simple result as 



follows, 



(i> s ,k=o\ exp ^ sA (A<f> at0 + A'n Sj0 )/Vl\ \i> a ,k=a) 
= exp - Si 



4pr]L y ' 4L. 



For charge sector, k = part of the initial density ma- 
trix is exp / ^/r t ^c.o) ' wnere -^co is the normal- 
ization TVco =Tr^cxp {^~P 2k\ n c oj ) ■ ^ nc evaluation of 
the charge sector yields 



-Tr < exp i 



cxp ( -p^Ln 2 



2K, 



c,0 



2 get 

2ttc c ^/3L 



Collecting the results above, we conclude 



exp 



4L 
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